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ABSTRACT 

While it has long been known that a large number of short-lived transient spirals can cause stellar 
migration, here we report that another mechanism is also effective at mixing disks of barred galaxies. 
The resonance overlap of the bar and spiral structure induces a nonlinear response leading to a strong 
redistribution of angular momentum in the disk. We find that, depending on the amplitudes of 
the perturbers, the changes in angular momentum, AL, could occur up to an order of magnitude 
faster than in the case of recurrent spirals. The signature of this mechanism is a bimodality in AL 
with maxima near the bar's corotation and its outer Lindblad resonance; this is independent of the 
properties of the spiral structure. For parameters consistent with the Milky Way the disk mixes in 
about 3 Gyr and the stellar velocity dispersion increases with time in a manner roughly consistent 
with observations. This new mechanism could account for both the observed age- velocity relation and 
the absence of age-metallicity relation in the solar neighborhood. Spiral-bar interaction could also 
explain observations showing that strongly barred galaxies have weaker metallicity gradients than 
weakly barred or non-barred galaxies. 

Subject headings: ISM: abundances — Galaxy: abundances — galaxies: ISM — Galaxy: evolution — 
galaxies: kinematics and dynamics — galaxies: structure 



1. INTRODUCTION 

The metallicity of the interstellar medium (ISM) in 
a galaxy is expected to increase with time due to the 
progressive enrichment through stellar feedback. Con- 
sequently, one would expect that younger stars formed 
at the same galactic radius would have higher metallici- 
ties. In addition, the ISM metallicity has been well estab- 
lished to decrease with increasing galactic radius, r, due 
to a more efficient star formation and enrichment of the 
ISM i n the central regions of g alaxies ([Wilson and Roodl 
Il994t iDaflon and Cunhall200l . This in turn would re- 
sult in coeval stars being progressively more metal-poor 
as r increases. 

However, the metallicities of stars in the solar neigh- 
borhood (SN) have b een observed to lack th e expected 
corre lation with age (|Edvardsson et all Il993t IHavwoodl 
2008). Re-examining the age-m etallicity distribution o f 
the Geneva-Cop enhagen survey ([Nordstrom et al.ll2004l ). 
IHavwoodl (|2008T) found that the data are only consis- 
tent with the age-metallicity relation (AMR) for stars 
younger than 3 Gyr. For other samples the width of the 
metallicity distribution was found to be larger than what 
could be accounted for by measurement errors and/or 
inhomogeneous ISM. Because of the age-velocity rela- 
tion (AVR), where older stellar samples are observed to 
have higher velocity d ispersion than younger ones (e.g. 
iHolmberg et al.l l2009f ) , stars from the inner and outer 
Milky Way (MW) disk can enter the SN and thus blur the 
expected AMR. However, this can only account for ab out 
50 % of the observed scatter (|Nordstrom et al.|[2003 ). 

A possible explanation for this discrepancy is the 
radial migration of stars from their birth gu iding 
radii ([Havwoodl 120081: ISchonrich and Binnevll2009D . In- 
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deed, there is now considerable evidence that the non- 
axisymmetry of the Galactic potential can cause signifi- 
cant perturbations in the motion of both stars and gas: 
observational evidence for this comes from (i) the ob- 
served non-circular m otions of gas flows in the inner MW 
([Bissantz et al.1 120031 ). (ii) the la rge non-axisym metric 
motions of star-forming regions (TXu et all [20061) . (iii) 
the moving groups in the SN ([Denned 11991" contain- 
ing stars of very different ages ([Famaev et al.l 120051 : 
lAntoia et al.l [2008). suggesting that their clumping in 
velocity space is most likely not due to irregular star 
formation , but ra t her to dynamical perturbati ons from 
the bar (iDehnenl l2000l: iMinchev et al.l l2010aD and/or 
the sp irals (e.g. jQuillen and Minchevil2005llAntoia et al.l 
20091), or an orbiting satellite galaxy ([Quillen et al.ll2009l : 
Minchev et aTll2009[ i. 

A numerically well-studied effect of this non- 
axisymmetry of the Galactic potential is the radial mi- 
gration of stars due t o resonant scattering by tra nsient 
spiral struc ture (SS) (ISellwood and Binnevl l2002t here- 
after SB02; IRoskar et al.ll2008D . which could explain the 
aforementioned absence of AMR in the SN. N-body sim- 
ulations based on cosmological initial conditions (e.g., 
iHeller et al.l [2001 would a priori be the most natural 
framework in which to study radial mixing in galactic 
disks. However, these simulations cannot at present pro- 
duce large enough thin galactic disks due to the trans- 
fer of angular momentum from the disk to the live 
dark matter halo. The resulting galaxies are typically 
thick and present large stellar velocity dispersions, so 
that they are closer to SO galaxi es than to MW-likc 
ones (lSanchez- Blazqucz ~et al.ll2009l ). On the other hand, 
IRoskar et al.l ([20081 ) used an N-body hydrodynamical 
model in which the initial conditions were designed to 
mimic the stage following a last major merger when the 
thin disk formation is supposed to start. Both SB02 and 
IRoskar et all (|2008h compare the effects of radial migra- 
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tion in the MW disk, for which the existence of a bar is 
firmly established, to a simulated galaxy lacking a central 
bar. Unless the effect of the bar is unimportant to the 
mixing process, their results would only poorly approx- 
imate the radial migration in our Galaxy. In this paper 
we check whether the effect of galactic bars can indeed 
be neglected in the context of radial migration in galactic 
disks, as assumed in a number of previous studies. For 
this purpose we study the combined effect of SS and a 
central bar (and more generally of multiple patterns) on 
the radial mixing of stars in a stellar disk. 

There is observational evidence for multiple patterns 
in external galaxi es, such as asymmetries in the sp i- 
ral structure (e.g. iHenrv et al.1l2003t iMeidt et al.ll2008D . 
By expanding galaxy images in Fourier components, 
lElmegreen et al.f (|T992) noted that many galaxies ex- 
hibit hidden three-armed components and suggested that 
multiple spiral density waves can propagate simulta- 
neously in galaxy discs. In addition, by m eans of 
the recently developed (jMerrifield et al.1 l2006f) Radia l 
Tremaine- Weinberg (TWR) method. IMeidt et al.l (|20Q9h 
analyzed the high-quality HI and CO data cubes avail- 
able for four spiral galaxies. The authors found direct 
evidence for the presence of resonant coupling of multi- 
ple patterns, including spiral-spiral and spiral-bar com- 
ponents. 

The results of N-body simulations provide additional 
motivation fo r studying the d yna mical effects of multi- 
le pa tterns. ISellwoodT(jl985| ) and ISellwood and Sparkd 
presented N-body experiments in which a bar co- 
exists with a spiral pattern movi ng a t a much lower angu- 
lar velocity. iTagger et all (|1987l) and lSygnet et al.l (|1988Q 
explained this as a non-linear mode coupling between 
a bar and SS. Sel lwood's simulations , and the theoreti- 
cal ex planation of Tag ger et al.l (|1987l) andlSygnet et al.1 
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w ere later confirmed by iMasset and Taggerl 
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How do multiple pattern s affect the dynamics of galac- 
tic disks? iQuillenl (|20Q3f ) considered the dynamics of 
stars that are affected by perturbations from both SS and 
a central bar by constructing a one-dimensional Hamilto- 
nian model for the strongest resonances in the epicyclic 
action-angle variables. Quillen pointed out that when 
two perturbers with different pattern speeds are present 
in the disc, the stellar dynamics can be stochastic, par- 
ticularly near resonances associated with one of the pat- 
terns. Simil ar findings were presented more recently by 
Ualalil ()2008( 1 . All these resu lts are not surpr ising since it 
has already been shown by lChirikovl (|19T9T ) that in the 
case of resonance overlap the last KAM surface between 
the two resonances is destroyed, resulting in stochastic 
behavior. We therefore expect that resonance overlap 
could give rise to both velocity dispersion increase and 
radial migration. 

In this paper we study the combined effect of a cen- 
tral bar and SS on the dynamics of a galactic disk. We 
concentrate on a steady-state SS in order to assess the 
pure effect of non-linear coupling between the two per- 
turbers, noting that while there is strong circumstan- 
tial evidence that SS is mostly not-stationary in un- 
barred galaxies, the situation is less clear in barred ones 
(|Binnev and Tremaind 12008'). In any case, the following 
model can also serve as a proxy for studying the effect of 
transients, and is complementary to the transient spiral 
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Fig. 1. — Resonances in a galactic disk for nearly circular or- 
bits and a flat rotation curve. Corotation occurs along the dotted 
black curve and is given by Hj, s = H(r), where ilf, s is the bar or 
spiral pattern speed, O(r) = vo/r is the local circular frequency, 
and vo is the constant circular velocity. The 2:1 outer and inner 
Lindblad resonances (OLR and ILR) occur along the solid black 
curves, computed as Q;, s = £l{r) ± ft/2, where re is the local ra- 
dial epicyclic frequency. The outer and inner 4:1 LRs occur along 
the dashed black curves, given by Jlj, s = H(r) ± re/4. The red 
horizontal line indicates the bar pattern speed used in this paper, 
£7;, = 55.5 km/s/kpc, and the blue horizontal line shows a spiral 
pattern speed of Q s = 24 km/s/kpc, such as the one used in the 
simulations shown in Figs. [3] and [4] The vertical red and blue lines 
give the radial positions of each resonance for the bar and spiral 
structure, respectively (solid lines: 2:1; dashed lines: 4:1, dotted 
lines: CR). 

mechanism developed by SB02. 

2. RESONANCES IN GALACTIC DISKS 

Galactic disks rotate differentially with nearly flat ro- 
tation curves, i.e., constant circular velocity as a function 
of galactic radius. In contrast, density waves, such as a 
central bar or SS, rotate as solid bodies. Therefore stars 
at different radii would experience different forcing due 
to the non-axisymmetric structure. Of particular inter- 
est to us are locations in the disk where the stars are 
in resonance with the perturber. The corotation res- 
onance (CR), where stars move with the pattern, occur 
when the angular rotation rate of stars equals that of the 
perturber. The Lindblad resonances (LRs) occur when 
the frequency at which a star feels the force due to a 
perturber coincides with the star's epicyclic frequency, 
k. As one moves inward or outward from the corotation 
circle, the relative frequency at which a star encounters 
the perturber increases. There are two values of r for 
which this frequency is the same as the radial epicyclic 
frequency. This is where the inner and outer Lindblad 
resonances (ILR and OLR) are located. Quantitatively, 
LRs occur when Qb.s = O ± n/m, where m is the mul- 
tiplicity of the pattern. The negative sign corresponds 
to the ILR and the positive to the OLR. While Bcrtil 
Lindblad defined these for the case of an m=2 pattern 
(thus strictly speaking the ILR/OLR are the 2:1 reso- 
nances), for an m=4 pattern the ILR/OLR must be the 
4:1 resonances. 

Since second order resonances, i.e., 4:1 for a 2-armed 
spiral or bar, or 2:1 for a 4-armed spiral, can also be quite 
important (as will be shown later) , we need to have a con- 
venient way to refer to them. It is somewhat confusing 
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and unclear how the 4:1 resonances are referred to in the 
literature. The inner 4:1 resonance for an m=2 pattern 
is known as the Ultra-harmonic resonance (UHR) . Some 
also describe the inner and outer 4:1 resonances as the 
IUHR and the OUHR, others as the inner and outer m=4 
resonance. If the pattern multiplicity is m=4, then these 
become the ILR and OLR. To our knowledge, there is no 
terminology for the 2:1 resonances of an m=4 pattern. 
To simplify things, in this paper we generalize the stan- 
dard notation of Lindblad resonances to allow for the 
existence of both 2:1 and 4:1 resonances, regardless of 
the multiplicity of the pattern. In other words, we will 
refer to the 2:1 ILR/OLR and the 4:1 ILR/OLR for both 
2-armed (or bar) and 4-armed SS. Naturally, other res- 
onances can also be described in this manner, i.e., 3:1, 
5:1, 6:1 ILR/OLR, etc. 

Fig. Q] illustrates the relationship between the pattern 
angular velocity and the radii at which resonances occur 
for a flat rotation curve and nearly circular orbits. Note 
that for two (or more) nonaxisymmetric patterns (such 
as bar + SS or SS + SS) moving at different angular 
velocities there will always exist regions in the disk where 
resonances overlap. As we will see later, this turns out 
to be crucial for the dynamics of barred galaxies. 

3. NUMERICAL PROCEDURE AND CHOICE OF 
PARAMETERS 

We follow the motion of test particles in an initially 
axisymmetric, 2D stella r disk with a flat rota t ion cu rve 
as described in detail in lMinchev and Quillenl (|2007l ). in 
which we grow non- axisymmetric perturbations as de- 
scribed below. Both the initial disk stellar density, S(r), 
and the radial velocity dispersion, ay(r), decrease ex- 
ponentially with radius, where we initially start with 
0>(?*d) = 5 km/s at the disk scale-length, = 3 kpc for 
a MW-type galaxy. To simulate the effect of a central 
bar and spiral density waves we impose p erturbations on 
the in itially axisymmetric disk as d o ne by Minchcv et al.l 
(|2007T) and iMinchev and Quillenl (|2007ft . respectively. 
Both bar and spiral perturbations are grown simultane- 
ously in 0.4 Gyr by increasing their amplitudes form to 
the maximum values. To ensure a smooth transition to 
the perturbed state we use the function defined in eq. 4 
from iDehnenl (|2000l ). To be able to relate parameters to 
the MW we work in units of the solar distance and cir- 
cular velocity with corresponding values ro = 8 kpc and 
v c = 240 km/s; t his gives flp = 30 km/ s/k pe. 

While in IMinchev and Quillenl (|2007l ) and 
IMinchev et al.l (|2007|) we required the integration 
of a large number of particles in order to be able to 
resolve a "solar neighbourhood", here we are interested 
in the global behavior of a galactic disk. Therefore we 
have found it sufficient to integrate 10 5 particles for 
each simulation, distributed between an inner and outer 
radii, Ti n = 2 and r out = 15 kpc. We have verified 
that increasing the resolution to 10 7 particles does not 
change our results. We divide each sample into 50 radial 
bins and calculate the azimuthally averaged changes in 
angular momentum in the disk in order to estimate the 
induced radial migration. The time evolution of each 
simulation is followed for 3 Gyr. We distinguish between 
the SS and the bar by indicating the OLR as OLR s and 
OLRb, respectively, and similarly for the case of the 
ILR. 



3.1. Perturbation from a central bar 

The bar strength, e&, for the case of the MW, is defined 
as the ratio of the forces due to the bar's quadrupole and 
the axisymmetric background at the Sun's ra dius, ro, 
on th e bar's major axis, e& = 3(At/vl)(rb/r ) 3 (|Dehnenl 
2000]), where v c is the asymptotic circular velocity. The 
bar potential is then given by 



$b = A h (e h ) cos[2{(f> - n h t)} x 




, r > r h 
, r < r h 



(1) 

where r& = 0.8r cr is the bar half-length in terms of 
its CR, r cr ; in physical units rf, = 3.44 kpc and r cr — 
4.3 kpc. 

Previous work ha s used as a measure of bar strength 
the parameter Qt (|Combes and Sanders! fl98lT ). This is 
the ratio of the maximum tangential force to the az- 
imuthally averaged radial force at a given radius. From 
cq.Q]this definition yields Qt = 2Ab/v^.. We examine bar 
amplitudes in the range 0.1 < Qt < 0.4 as expected from 
observa tions of various galaxies an d from N-body simu- 
lations (jCombes and Sa nders 1981). This corresponds to 
0.013 < | et, | < 0.05 in our units. 

The bar pattern speed is chosen to be fif, = 55.5 
km/s/kpc, which for the case of our Galaxy can be writ- 
ten as fib — 1.85f2o- This value is consistent with obser- 
vations of external galaxies where it is foun d that bars are 
dyna mically fast, i.e., r cr /n> = 1.2 ± 0.2 (jAguerri et al.1 
l2003h . and is also in accordance w ith recent bar pattern 
speed determinations for the MW ( Minchev et al.ll2010al 
I2007t iBissantz et al.ll200llDehnenll2000l ). 

3.2. Perturbation from spiral structure 
The spiral potential is given by 

r 

& a (r,<f>,t) = e s cos[Q;ln m(<f> - O s t) , (2) 

ro 

where e s is the SS strength, related to the amplitude of 
the mass surface density of spirals, S s , as 

e s « -2irG^ s r /{av 2 c ), (3) 

as shown in iBinnev and Tremaind ([2008). The param- 
eter a is related to the pitch angle of the spirals, p, by 
a = rncot(p). The azimuthal wavenumber m is an inte- 
ger corresponding to the number of arms. We consider 
both 2-armed and 4-armed SS with a = —4 and -8, re- 
spectively, where the negative sign corresponds to trail- 
ing spirals. lElmegreenl (|1998l) found that grand-design 
spirals have arm-interarm contrasts of 1.5-6, correspond- 
ing to a fractional am plitude of 0.2 < S „ /£ < 0.7, which 
is in agreement with iRix an d Zaritsky (1995) who esti- 
mated 0.15 < (£ s /£) < 0.6. Since we would like to ex- 
amine this range of spiral amplitudes we need to relate 
the relative overdensity to the relative potential. 

For a maximum exponential disk the peak circular 
speed in the disk has the value v c ~ 0.622 \JGMdJr~d 
at r ~ 2.157^, where is the disk mass inclosed 
by r^. The surface density of the disk at radius r is 
£(r) = (Af ( i/27rr^)e _r / rd . Eliminating Md from these 
expressions, we find at r , v% ~ 0.39 x 2-KGrdY>oe ra l rd . 
Substituting this expression for v c in eq. [3l the relation 
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Fig. 2. — Changes in angular momentum, AL, as a function of the initial angular momentum, Lq. Both AL and Lq are in units the 
circular velocity, v c ; Lq is thus equivalent to the galactocentric distance, r, in kpc. Bar and SS amplitudes are indicated in each panel. The 
first two panels form left to right present simulations with a bar of an intermediate strength (Qt = 0.25) and a strong one (Qt = 0-4), 
respectively. The third and fourth panels show a 4-armed SS with relative overdensity S s /So = 0.2,0.35. The rightmost panel presents 
a simulation with a 2-armcd SS with an amplitude giving rise to £ s /Eo = 0.25. The dotted lines show the corotation radii. The 2:1 and 
4:1 LRs are indicated by the solid and dashed lines respectively (bar=red, spiral=blue). Note that, depending on the pattern speed, some 
resonances might not be present in the disk. |AL| increases significantly only near the corotation of each perturber. 
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Fig. 3. — Time evolution of the changes in angular momentum, AL, for a stellar disk perturbed by both a central bar and SS. The top 
row shows a bar and a 4-armed SS of intermediate strengths: e\y = —0.03 (Qt = 0.25) and e s = —0.02 (E s /Eq = 0.35 for m=4). For the 
bottom row the bar and the 2-armed SS are weaker: = —0.02 (Qt = 0.16) and e s = —0.03 (E 3 /Eo = 0.25 for m=2). The increase of 
AL with time indicates that stars are being placed on radii different than their birthplaces, i.e., radial migration takes place throughout 
the disk. 

between the relative potential and the relative overden- 
sity becomes 

ts ~ S r d 0.39a' 1 ' 

To correct for finite disk thickness we need to de- 
crease the righ t-hand side of this equation by the fac- 
tor e~ kz ° (e.g., I Vandervoortl 1 1 9 70l) . where k is the wave 
vector given by k rs a/ro. For a MW-type galaxy we 
can adopt a disk scale-height of zq ~ 300 pc, resulting 
in e~ kza » 0.86 and 0.74 for a 2- and 4-armed SS, re- 
spectively. Using eq. HI and accounting for finite disk 
thickness we estimate that the observed spiral arm am- 
plitudes quoted above would result in potential pertur- 
bation in the range 0.007 < |e s | < 0.031 for a 4-armed 
and 0.015 < |e s | < 0.072 for a 2-armed SS. 

4. SINGLE PERTURBATION 

We now estimate the changes in angular momentum, 
AL, in a stellar disk. To do this, we divide the disk into 
50 radial bins and compute AL = L\ — Lq, where Lq and 
L\ are the mean values of the initial (at t — 0) and final 
(at some later time, t\) angular momenta of particles in 
each bin. 



In Fig. [2] we show AL as a function of Lq for a simula- 
tion with a single perturber. Both axes are in units of the 
circular velocity, v c ; Lq is thus equivalent to the galacto- 
centric distance, r in kpc. The bar and SS amplitudes are 
indicated in each panel. The first two panels from left to 
right present simulations with a bar of an intermediate 
strength (Qt = 0.25) and a strong one (Qt = 0.4), re- 
spectively. The third and fourth panels show a 4-armed 
SS with relative overdensity £ s /£ = 0.2,0.35. The 
rightmost panel presents a simulation with a 2-armed SS 
with an amplitude giving rise to S s /So = 0.25. The dot- 
ted lines show the corotation radii. The 2:1 and 4:1 LRs 
are indicated by the solid and dashed vertical lines, re- 
spectively (bar=red, spiral=blue). Note that, depending 
on the pattern speed, certain resonances might not be 
present in the disk. The pattern speeds are given in units 
of the angular velocity at r = 8 kpc, O = 30 km/s/kpe. 

Significant angular momentum changes in Fig. [2] occur 
only near the CR of each perturber (dotted lines). No 
further increase in \AL\ is apparent with time once the 
perturber is fully grown, which is to be expected for a 
time-independent perturbation. We are now interested 
to see what changes occur when both the bar and SS act 
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n„=i.85, n s =o.80 

e t ,= T 0.01,e s =-0.01 
m s =:4 I I 



Q b =1.85, Q s =0.80 
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Q b =1.85, Q s =0.80 
~e h =-0.03,e =-0.01 



Q„=1.85, Q s =0.80 
e,=-0.Q4,e =-0.01 



0.85, n s 
-0.05j 



= 0.80 
= -0.01 




together. 



5. BAR + SPIRAL STRUCTURE 



Fig. 4. — The changes in angular momentum at t = 3 Gyr for simulations with different perturbation strengths, while all other parameters 
are kept fixed. From left to right the bar strength increases as indicated in each panel. The SS strength increases from top to bottom. 
Note the bimodality caused by the bar's corotation and its 2:1 OLR. 

SS have intermediate strengths: ej, = —0.03 (Qt = 0.25) 
and e s = -0.02 (E s /S = 0.35 for m=4). For the bottom 
row the bar and the 2-armed SS are weaker: e& = —0.02 
(Qt = 0.16) and e s = -0.03 (E s /S = 0.25 for m=2). 
The pattern speeds in both simulations are the same so 
that we can asses the effect of the different amplitudes. 

Examining the first row of Fig. we note that at 
the beginning of the simulation |AL| increases mainly 
at the CR of each perturber (dotted lines), which would 
be the case of linearly adding the individual effects of 
the bar and SS seen in Fig. [5] However, at later times 
large changes in angular momentum occur, suggesting 
that non-linear effects become important. Note that at 
t = 0.6 Gyr the effect near each CR is already stronger 
than that expected from each individual perturber (com- 
pare to Fig. [5]). Meanwhile the area between the CRs of 
the bar and SS fills in with time as a result of the proxim- 
ity of the 4:1 ILR S , the 4:1 OLR b and the 2:1 OLR b , thus 
extending the effect over the entire disc. At t = 1.8 Gyr 
(fourth panel of Fig. [3]) , the changes in angular momen- 
tum throughout the disk are similar in amplitude to those 
resulting from recurrent spirals (Fig. 11 in SB02). To 
reach this level of mixing, however, requires less than 
~ 2 Gyr compared to ~ 9 Gyr for the case of transients. 

Decreasing the bar strength to e& = —0.02 (Qt = 0.16) 
(bottom row of Fig. [3J and considering a slightly weaker 



As shown by Mi nchev and Qui llcn (2006J, when two 
perturbers moving at different pattern speeds are im- 
posed on a galactic disk, an increase in the random mo- 
tions of stars is expected, i.e., the disk heats. In addition, 
the authors found that in regions of resonance overlap 
stars drift radially with time from their birth radii. This 
suggests that such a mechanism could be responsible for 
radial migration in the disk. Here we quantify this idea 
more clearly for the case of bar + SS, by examining the 
effect of the simultaneous propagation of the two per- 
turbers. 

5.1. Time evolution of AL 

In Fig. [3]we show the time development of the changes 
of angular momentum in a stellar disk perturbed by both 
a central bar and SS. In each row we present a different 
combination of bar and SS strengths, as described bellow, 
while all other parameters are kept fixed. As in Fig. [21 
contour plots show the change in angular momentum, 
AL, as a function of the initial angular momentum Lq 
(or radius, r in kpc). Panels from left to right show the 
temporal evolution of the system for 0.3, 0.6, 1.2, 1.8 and 
2.5 Gyr. In the top row both the bar and the 4-armed 




Fig. 5. — Variation of mixing with pattern speed ratio. AL is shown for different SS pattern speeds, Q s , while keeping Qi, the same. 
This results in different spacing of the bar and SS resonances. From left to right the 4-armed SS moves at Q a = 0.7, 0.8, 0.9, 1.0, and 1.1 in 
units of f!rj = 30 km/s/kpc. All other simulation parameters are kept fixed: the time is t = 3 Gyr and the strengths of the perturbers are 
the same as in the top row of Fig. [3] As expected, the effect of resonance overlap on AL is strongest when the highest order resonances 
(2:1 OLRt, and 4:1 ILR S ) lie closest to each other (leftmost panel, Lq ~ 7.5). 



2-armed SS has a weaker effect on AL. However, changes 
in AL are still very strong, mixing the disk in less than 
3 Gyr. 

At this time we can already identify two main differ- 
ences between radial migration from transients and the 
mechanism described in this work: (1) The time-scale 
for disk mixing is much shorter in the case of resonance 
overlap and (2) resonance overlap causes a bimodal dis- 
tribution in AL with peaks near the bar's CR and its 
2:1 OLR, compared to the single peak seen in Fig. 11 
and 12 in SB02. Next we explore further the effect of 
perturbation strengths. 

5.2. Variation with perturbation strengths 

The strong increase of | AL| with time when both a bar 
and SS perturb a stellar disk suggests that this is a non- 
linear effect. We therefore expect strong variation with 
the strength of the perturbers. To investigate this fur- 
ther we now explore the effect on mixing when different 
combinations of e s and are considered. In Fig. [4] we 
plot the changes in angular momentum at t = 3 Gyr for 
simulations with different perturbation strengths, while 
all other parameters are kept fixed. From left to right the 
bar strength increases as indicated in each panel. The SS 
strength increases from top to bottom. 

Observing the top, leftmost panel of Fig. @] we see that 
when both perturbers are weak mixing is not very strong, 
although an increase in |AL| near the 2:1 OLRt, is still 
apparent. However, as the bar strength increases, while 
keeping the spiral weak, the resonance overlap in the in- 
ner disk has a surprisingly strong effect on AL. To realize 
this compare the top row of Fig. [4] to the corresponding 
bar strengths in Fig. [5J Similar behavior is observed in 
the case of a weak bar and a range of SS strengths (first 
column in Fig. |4j . Note that here in addition to the inner 
disk we also see a strong effect near the spirals' CR at 
10 kpc, which is well outside the 2:1 OLRt,. A possible 
explanation could be the proximity of the CR of SS to 
the bar's 1:1 asymme tric , 1:1 symm etric resonances, or 
x\{2) unstable orbits ()Dehnenll2000l ). 

5.3. Variation with pattern speed ratio 

In Fig. [5] we plot AL for different SS pattern speeds, 
fl s , while keeping flj, the same. This results in different 
spacing of the bar and SS resonances. From left to right 
the 4-armed SS moves at n s = 0.7,0.8,0.9,1.0, and 1.1 
in units of f2o = 30 km/s/kpc. All other simulation 
parameters are kept fixed: the time is t — 3 Gyr and 
the strengths of the perturbers are the same as in the 



top row of Fig. [3] As one would expect, the effect of 
resonance overlap on AL is strongest when the highest 
order resonances (2:1 OLRb and 4:1 ILR S ) lie very close 
to each other (leftmost panel, Lo ~ 7.5). For all pattern 
speed ratios the distribution of AL exhibits a bimodality 
near the bar's CR and near its 2:1 OLR. This signature 
can be used to identify the mechanism described here in 
fully self-consistent simulations, where the disk dynamics 
is much more complicated and may include additional 
mixing due to transient spirals and small mergers. 

6. IMPLICATIONS FOR THE MILKY WAY 

We showed in Sec.[5]that the simultaneous propagation 
of a central bar and SS in a stellar disk induces radial 
migration. Since this is a nonlinear effect, depending on 
the perturbers' amplitudes, strong mixing could occur 
in a short period of time. Similarly to the majority of 
disk galaxies, the MW has been well established to be a 
barred galaxy, thus we expect mixing due to resonance 
overlap to be at work here as well. While there is no 
doubt that the MW has been affected by this process, it 
may be difficult to estimate exactly how. This migration 
mechanism is a strong function of the strengths of the 
MW bar and SS. However, it may be incorrect to use 
the currently observed spiral and bar amplitudes since 
most likely they have not been the same throughout the 
lifetime of the Galaxy. In fact, we expect that spirals 
were stronger in the past, when the disk velocity disper- 
sion was lower, thus resulting in a more gravitationally 
unstable disk. In any case, given that the past history of 
the MW spiral and bar structure is largely unknown, we 
will make an estimate of how much mixing could have 
resulted by using the currently observed SS and bar pa- 
rameters. 

6.1. Choice of bar and spiral structure parameters 

For all simulations in this paper we have used a 
bar pattern speed of fib — 1.85Qq, where f^o is the 
angular velocity at the solar circle rrj. For ro = 
8 kpc and v c — 240 km/s, in physical units fib — 
55.5 km/s /kpc, consistent with recent pattern speed 
estimates (iMinchey et all [2010al 120071: iBissantz et al. 
l200aiDehnenll200Q[ ). iRodriguez -Fernandez and Combed 
(|2008h estimated the strength of the MW bar to lie in 
the range 0.25 < Qt < 0.4 by fitting bar models to near 
infrared observations (2MASS counts). In the units used 
in this paper this range corresponds to 0.03 < \eb\ < 0.05. 
The individual effect on AL of bars with amplitudes 
£h = —0.03 and e& = —0.05 can be seen in Fig. [2] 
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Fig. 6. — The distribution of the initial radii of stars, which have ended up in an annulus of 400 pc (in the green bin) centered at a given 
galactocentric distance, r. Simulation parameters are consistent with the Milky Way structure (see Sec.[6]l. From top to bottom different 
rows show the time evolution of each such annulus centered on r = 2, 4, 6, 8, and 10 kpc. The black and red histograms show stars on 
nearly circular orbits (ay ~ 5 km/s) and the total population, respectively. With this new mechanism, the solar neighborhood stars could 
have come from a wide range of Galactic radii, including from the outer Galaxy, and thus account for the flatness of the observed AMR. 
Note the different signature in the mixing of cold and hot orbits. 



The MW SS parameters are much more uncertain. It 
has been proposed to be a superposition of 2- and 4- 
armed patterns moving near corotati on with the Sun 
(iLepine et al.1 12001b iTorra et all l2000ft. On th e other 
hand. iMartos et all 120041 ) and iBissantz et al.l (|2003( ) 
matched the properties of the gas in the inner Galaxy 
with a 4-armed SS with a pattern speed of 20 km/s/kpc 
(or with the Sun near the 4:1 ILR). The MW SS 
at this stage in the Galaxy evolution appears to be 
quit e weak compared to observa tions of external galax- 
ies: iDrimmel and Spergell ((2001h derived S s /So = 0.15 
for the stellar component (based on K-band observa- 
tions). These authors, however, suggested that most 
likely this value is larger. There are also indications 



from the gas flow in the inner Galaxy that the MW has 
stronger spiral arms, where the amplitude in the mass 
density is larger by a factor of 1.5 than it s amplitude in 
the near-infrared luminosity den sity (see IBissantz et al.l 
120031 : iFamaev and Binnevl 12005) . Therefore, a possible 
range for the relative overdensity is 0.15 < S s /So < 0.23. 
These values, including the correction for finite disk 
thickness (Sec. I3.2[) . correspond to a potential pertur- 
bation of 0.015 < |e s | < 0.024 and 0.007 < |e s | < 0.01 for 
a 2- and 4-armed SS, respectively. 

From the discussion above, for parameters consistent 
with our Galaxy, the amount and rate of mixing due to 
spiral-bar resonance overlap is given by the second row of 
Fig. |U for a 4-armed SS, and by the bottom row for a 2- 
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Fig. 7. — Radial velocity dispersion as a function of time at the 
solar radius (solid line) for the simulation in Fig. \E\ The dashed 
line shows a power-law fit to the simulated data. The dotted line 
plots the observed age -velocity relation in the solar neighborhood 
l|Holmberg et alj|2009h . 

armed SS (excluding the leftmost panels in both cases). 
Note that the strength of the 2-armed SS in Fig. 2] is 
slightly stronger than the estimated range. 

6.2. Mixing at different Galactic radii 

To see in greater detail how stars migrate throughout 
the disk we now examine the effect on mixing at different 
galactic radii. To be consistent with MW parameters, we 
use the simulation for which we presented the changes in 
angular momentum in Fig. [4l second row, second col- 
umn. In Fig. [6] we show the distribution of birth radii 
(at t = 0) of stars which are "currently" (at the time 
indicated in each plot) in an annulus of 400 pc (in the 
green bin). Each row shows an annulus centered on a 
different galactic radius, r: from top to bottom r = 2, 
4, 6, 8, and 10 kpc. To show the effect on stars placed 
on nearly circular orbits, the black histograms plot only 
stars with a cut in the eccentricity, e < 12. We estimate 
e as e ~ \/(u 2 + 2u 2 )/2, approp riate for a flat rotation 
curve ijArifvanto and Fuchsll2006l ). and u, v are the radial 
and tangential velocity components of stars with the cir- 
cular velocity subtracted. The radial velocity dispersion 
resulting for this sample is ~ 5 km/s, insuring that stel- 
lar orbits are almost circular. The red histograms show 
the full sample. Stars on cold orbits, shown by the (nor- 
malized) black histograms, constitute about 10% of the 
total population. 

Firstly, note that the distributions peak at radii near 
regions of resonance overlap, as expected. While in the 
inner disk (r = 2 and 4 kpc) there are strong weights at 
larger radii, the situation is reversed for r > 6 kpc. At 
the solar distance, r = 8 kpc (fourth row in Fig. [6]), a 
large population of stars from the bar's CR (r ~ 4.5 kpc) 
enters the solar circle. Naturally, when the total popula- 
tion is considered (red histogram) more stars enter from 
the inner disk due to the exponential profiles of both the 
density and the radial velocity dispersion. Observe also 
the overall shift toward larger radii of the cold population 
compared to the total sample. 

While the pattern speed of the MW SS is not well 
constrained, the bar is well established to have its 2:1 
OLR just inside the solar circle. Since the bar's CR and 
OLR are both strong and roughly invariable (only small 
changes in the bar parameters are expected in ~ 3 Gyr) , 
we can expect to find a bimodal distribution of AL with 
peaks at these bar resonances regardless of the SS pattern 



speed. That is exactly what is seen in Fig. [5] What if the 
MW spirals were transient? In that case the distribution 
of the changes in angular momentum can be approxi- 
mated by the linear combination of a range of simulations 
with different SS pattern speeds (imagine adding up all 
panels in Fig. |5|) . This would again result in a bimodal 
distribution. Therefore, although the past history of the 
MW structure is not known and a large range of pos- 
sible perturber strengths, pattern speeds and longevity 
exists, it would be reasonable to expect to always find a 
bimodality in the distribution of AL. 

Consequently, the large fraction of stars coming from 
the bar's CR which enters the solar circle (fourth row 
of Fig. may be one of the observational signatures of 
the spiral-bar resonance overlap mechanism. This is ex- 
pected to affect the local metallicity distribution of stars. 
We also note that a large fraction of stars comes from the 
outer disk, roughly 40% of all stars for this particular 
pattern speed ratio. This is consistent with the metal- 
poor stars from the outer Galaxy scattered int o the SN, 
as re cently inferred from observational data (|Havwoodl 
2008), and a metal-rich population coming from the inner 
Galaxy. With this new mechanism, the SN stars could 
have come from a wide range of Galactic disk radii, in- 
cluding from the outer Galaxy, and thus can account for 
the flatness of the observed AMR. 

On the other hand, away from the solar circle we expect 
to find overdensities of stars originating at the bar's 2:1 
OLR. We clearly see this in the initial distribution of 
stars ending up in the annuli centered on r = 2, 4, 6 and 
10 kpc (Fig.[S]). This may provide another observable for 
identifying the currently discussed mechanism. 

However, the most promising technique to put con- 
straint s on this mechanism in the MW is " chemical tag- 
ging" (|Freeman and Bland-Hawthornll2002f h The idea is 
to use the detailed chemical abundance patterns of indi- 
vidual stars in our Galaxy to associate them with com- 
mon ancient star-forming aggregates with similar abun- 
dance patterns. This technique will be very important 
to characterize the star formation and chemical histo- 
ries of the MW. Chemical tagging will become possi- 
ble with the HER MES survey on the Anglo-Aus t ralian 
Telescope (AAT; [Freeman and Bland-Hawthorn! 120081 : 
iFreeman et all l2010f ) and t he APOGEE survey at th e 
Apache Point Observatory (|Allende Prieto et al1l2008D . 
Therefore, in the next couple of years we may be in a 
position to directly measure the spread of open clusters 
across the Galaxy as a function of age, thus providing 
direct constraints on the amount o f mixing in the Milky 
Way (|Bland-Hawthorn et alJl201Ch . 

6.3. The effect on disk heating 

Now we look at the effect on the increase of ran- 
dom motion of stars with time. The observed cor- 
relation between the ages and velocity dispersions 
of SN st ars has been a subject o f study sinc e the 
work of ISpitzer and Schwarzschildl (|1951l |1953|) (see 
Minc hev and Quillenl 120061 for a detailed introduction). 
Deferent mechanisms have been proposed to explain 
the age-velocity relation (AVR) (or stellar disk heat- 
ing) , including scattering by giant mole cular clouds 
(GMC) (e.g., iMihalas and Binnevl 119811) and tran- 
sient spirals (e.g., Barbanis and Woltierl Il967). In 



iMinchev and Quillen ( 20061 ) we showed that the reso- 
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nance overlap of multiple spiral patterns provides yet 
another mechanism for disk heating. We therefore ex- 
pect that an overlap of bar and SS resonances would 
give a similar result. To check if this is indeed the case, 
in Fig. [7] we plot the radial velocity dispersion, ay, (solid 
line) as a function of time for the same simulation used 
to produce Fig. [5] The dashed line shows a power-law 
fit to our simulated data: oy ~ t 29 . To compare to 
observations we also show the AVR in the SN (dotted 
line), as estimated recently from analysis of the G eneva- 
Copenhagen Survey data (jHolmberg et al.l [20091) . Note 
that the observed AVR goes like oy ~ i 39 . Moreover, in 
our simulation the disk heats ~ 20% more than observed. 
There could be several possibilities that could account for 
this discrepancy: (1) our simple model does not include a 
gaseous disk component nor star formation, which would 
naturally lower the velocity dispersion; (2) while here we 
considered 2D disks, a significant fraction of the heat- 
ing in the galactic disk plane would be transferred to the 
vertical direction, where GMC serve as s cattering agents 
(jJenkins and Binnevlll99 0: Jenkins 1992); and (3) in this 
work we used constant parameters for both the MW bar 
and SS, but it is reasonable to expect more realistic sce- 
narios, where throughout the lifetime of the Galaxy there 
could have been variations of the perturbers' amplitudes 
and pattern speeds, bar size, spir al pitch angle, as well as 
bar d estruction and reformation (jBournaud and Combed 
I2002T) . 

7. DISCUSSION AND CONCLUSIONS 

In this work we have examined the combined effect of 
a central bar and spiral structure on the radial migration 
in galactic disks. We have found that: 

(i) This new mechanism for galactic disk mixing could 
be up to an order of magnitude more effective than the 
transient spirals one (see Sec. I5.2|) . We note that the 
effect is non-linear, strongly dependent on the strengths 
of the perturbers. The non-linearity of the mechanism 
is confirmed by the fact that an increase in the ampli- 
tudes yields a much larger increase in the angular mo- 
mentum changes in a simulation including both bar and 
spiral structure, compared to the case of a single per- 
turber. In other words, the individual effects do not add 
up linearly. Another indication of non-linearity is that 
the effect increases with time only in the case of the two 
perturbers propagating simultaneously. 

(ii) The signature of this mechanism is a bimodality in 
the changes of angular momentum in the disk with max- 
ima near the bar's corotation and its outer Lindblad res- 
onance (see Figs. [5] and [6]) . This is true regardless of the 
spiral pattern speed and can be used in identifying the 
spiral-bar resonance overlap migration in N-body simu- 
lations, where the disk dynamics is much more compli- 
cated. 

(iii) For bar and spiral parameters consistent with Milky 
Way observations we find that it takes ~ 3 Gyr to achieve 
the mixing for which transients require 9 Gyr. In addi- 
tion to radial mixing, spiral-bar coupling can account 
for the age- velocity relation (AVR) observed in the solar 

2 Note that resonance overlap, a nd thus migration, wi ll also re- 
sult from multiple spiral structure l lMinchev et aljl2010bf ). 



neighborhood (Fig. Sec. I6.3[) . Note however that the 
estimates for the Milky Way spiral and bar strengths de- 
rived from observations now may simply be irrelevant for 
our understanding of the current state of mixing in the 
Galactic disk (see Sec. [6]). 

(iv) This mechanism could explain observations showing 
that strongly barred galaxies have weaker metallicity gra- 
dients than weakly barred or non- barred galaxies (e.g., 
I Zaritskv et all 119941 : iMartin and Royl 11994 [Perez et alJ 
2009). 

The effect of resonance overlap we have described in 
this paper appears to be quite strong. In order for it to 
work, a galactic disk must contain a central bar 2 . Unlike 
in our test-particle simulations, where the disk responds 
solely to the bar potential, as seen in N-body simula- 
tions and deduced from observations, a bar is expecte d 
to always drive spiral structure (e.g., ISalo et all 12010). 
Thus, any barred galaxy would be affected strongly by 
this mechanism. As the disk heats and SS diminishes 
in strength, the effect on AL will decrease. It is there- 
fore likely that resonance overlap migration will be most 
vigorous at certain periods during a galaxy lifetime. 

In this work we have only considered steady-state spi- 
rals in order to assess the pure effect of the coupling 
between bar and spiral structure. Considering tran- 
sient spirals would increase the efficiency, since the SB02 
mechanism would come into play as well. Therefore, a 
combination of a central bar and short-lived transients is 
expected to be extremely efficient at mixing barred spiral 
galaxies. 

As discussed in Sec. O the most promising tech- 
nique to put constraints on this mechanism in the MW 

is " ch emical tagging" (|Fr ccman and B land- Hawthorn! 

120021 iBland-Hawthorn et all l2010f > which will be- 
come possible w ith the forthcoming spectroscopic 
survey HERMES rtFreeman and Bland- Hawthorn! [20081 : 
iFreeman et~aT1l2010[ ). More ongoing and planned large 
Galactic surveys, such as RAVE, SEGUE, SIM Lite, 
GYES, LAMOST and GA IA, can search for signatur es 
of the mechanism (see also Minchev and Guillen 2008) . 

In a follow-up paper (jMinchev et al.ll2010 b) we study 
the effect of resonance overlap in fully self-consistent, 
TreeSPH/N-body simulations including a gaseous com- 
ponent and star formation. The migration mechanism 
described in this paper can provide an explanation for 
the observed flatness and spread in the age-metallicity 
relation. Our results, however, need to be incorporated 
into a chemo-dynamical galactic model to properly as- 
sess the effect on the AMR, and possibly also on the 
formation of the thic k disk , similarly to the work by 
ISchonrich and Binnevl (|2009l ). 
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